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The three-dimensional evolution of a plane wake 

By H. Maekawa, 1 R. D. Moser 2 AND N. N. Mansour 2 
1. Introduction 

In the past three decades, linear stability analysis has led to a comprehensive 
understanding of the linear stages of transition in plane wakes. Our understanding 
of the nonlinear and turbulent stages is less developed. Nonlinear theory developed 
by Papageorgiou k Smith (1988) was used to study the long- wavelength regime in 
wakes. The nonlinear and turbulent stages have been investigated experimentally, 
and few numerical studies have examined the early nonlinear stages of forced wakes. 

Experimental studies of the wake have been carried out extensively by Sato k 
Kuriki (1961), Sato k Onda (1970), Sato k Saito (1975, 1978), and Gharib k 
Williams-Stuber (1989). Sato k Onda (1970) reported that the wake responds 
differently under different kinds of forced disturbances. Gharib k Williams-Stuber 
(1989) investigated the enhancement and cancellation of perturbations in a plane 
wake by the strip heater technique. They showed that the mean velocity profile 
adjusts itself to become more receptive to the forced frequency. 

In their numerical study, Maekawa, Mansour k Buell (1992, hereafter referred 
to as MMB), simulated spatially developing plane wakes with inlet forcing using a 
single frequency (fundamental mode), a fundamental mode with its subharmonics 
and random perturbations. They found that phase jitter around the fundamental 
frequency plays a critical role in generating vortices of random shape and spac- 
ing. They also observed pairings and doublets in the wake. These observations 
were similar to those reported by Comte, Lesieur k Chollet (1989) and by Aref 
k Siggia (1981) for temporally developing layers. Furthermore, MMB found that 
infinitesimal random disturbances generate clear alternate-signed vortex streets due 
to selective amplification of the fundamental mode. 

The formation of three-dimensional structures and turbulent breakdown are not 
as well understood. Chen, Cantwell k Mansour (1990) showed that the phase 
between a two-dimensional fundamental mode and oblique modes controls the for- 
mation of three-dimensional structures in plane wakes. Depending on the phase, 
loop-like structures were observed. These structures were similar to the flow visu- 
alization experiments by Breidenthal (1980), Meiburg k Lasheras (1988), Antonia, 
Brown k Bisset (1987), and Perry, Lim k Chong(1980). However, a detailed under- 
standing of the underlying mechanism that initiates three-dimensional breakdown 
in wakes and the resultant structure is still missing. 

In what follows, the evolution of three dimensional disturbances in a incompress- 
ible wake is investigated using direct numerical simulations. The instantaneous 
three-dimensional structures and corresponding statistics are presented. 
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2. Preliminaries 

The time-developing mixing layer code used by Rogers Sc Moser (1992) was mod- 
ified to accommodate wake profiles. The cases described in this study are initialized 
with controlled initial conditions. We consider temporally evolving layers where the 
flow is taken to be periodic in x- and z- (streamwise and spanwise) directions. The 
cross-stream direction, y, is taken to be doubly infinite with boundary conditions 
applied at x — ±oo. The length of the computational domain in the streamwise 
and spanwise directions, L z and L : , is set to be integer multiples of the perturba- 
tion wavelength, Thus, in general, L x = n\ x and L z = m\ z . n and m are one or 
two for all cases reported here. The Reynolds number defined by the half-width 
and the deficit velocity of the wake is 300. As was done in MMB, the flow struc- 
ture and various flow statistics were examined for a range of Reynolds numbers, 
100 < Re < 600. The observed variations with Re indicate that Re = 300 is large 
enough to eliminate low-Reynolds number effects. 


2.1 The governing equations and numerical considerations 
The simulations reported here were performed by solving the vorticity equation 
derived from the incompressible Navier-Stokes equations: 

W + v * ( " xt,) = K vV (1) 

where Re o — U^bi^/v (v is the kinematic viscosity) is the Reynolds number, 2>i / 2 
is the initial half-width of the wake, and U 0 is the centerline deficit velocity. All 
quantities are non-dimensionalized by the appropriate characteristics scale, Uo and 

&i /2 • 

The spectral Galerkin method described by Spalart, Moser Sc Rogers (1991) 
was used for spatial discretization. The equations were advanced in time using a 
compact third-order Runge-Kutta scheme. A typical three-dimensional simulation 
(n = l,m = 1) required 48 x 64 x 48 Fourier/ Jacobi modes. 

2.2 Specification of initial conditions 

The initial mean velocity used for all the simulations reported here is of the form 

rr = cxp(-ln(2)(^) ) (2) 


which has a half width of 6 1/2 * In addition to the mean flow, simple perturbations 
are included in the initial conditions. 

To specify the initial conditions and facilitate discussion through this paper, it is 
necessary to refer to specific wavenumbers; they will be referred to in ordered pairs 

/ f k x X x k z X z ^ 


(3) 


where X x and X t are the x and z wavelengths from linear stability theory, k x and 
k z are corresponding wavenumbers. The fundamental modes discussed here have 
a and/or /3 equal to one. The amplitude A a p of a Fourier mode (a,/?) is defined 
to be the square root of the integrated (in y) energy in that mode, as described in 
Rogers Sz Moser (1991). 
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FIGURE 1. Spanwise vorticity structure of the two-dimensional wake (Case 2D-1) 
at the first maximum of A\ 0 (t = 48). Contour increments are 0.05 and tic marks 
are at b\j 2 intervals. 



FIGURE 2. Time development of the two-dimensional fundamental mode amplitude 
for Case 2D-2 (high initial amplitude) and 2D-1 (low amplitude). 

3. Accomplishments 

3 A Two-dimensional simulations 

To facilitate discussion of the three-dimensional evolution of the wake, we start by 
briefly reexamining the behavior of two-dimensional wakes. The spatially developing 
simulations of MMB showed that the growth of the fundamental mode is responsible 
for the formation of the vortex street in wakes. We consider three different initial 
conditions for the current two-dimensional simulations. Case 2D-1 starts with a 
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single fundamental mode with small amplitude (A\o = 0.01), Case 2D-2 starts 
with a single fundamental mode with large amplitude (A\o = 0.5), and Case 2D-3 
starts with a fundamental mode and its first subharmonic with large amplitudes 
(A \ o = Aio = 0.5). 

As in the spatially developing case, the amplitude of the fundamental grows 
exponentially, then saturates. Contours of spanwise vorticity at the time when the 
amplitude of the fundamental disturbance, reaches its maximum are shown 

in Figure 1 for Case 2D-1. The vortices evident in this figure have an oblong or 
elliptical shape, with the major axis roughly vertical. Successive snapshots show 
that the structures rotate by 90° so that the major axis is horizontal, consistent 
with the observations of MMB. If we define the rollup time t r to be at the first 
maximum of the fundamental disturbance amplitude, we find for this case, t r = 48. 

When a larger initial amplitude is used (Case 2D-2), the fundamental mode sat- 
urates sooner (* r = 18). After saturation, the mode loses energy at a much faster 
rate than in Case 2D-1, then starts growing again, and a second maximum appears 
(see Figure 2). The orientation of the major axis of the vortical structure at the first 
maximum of A\ 0 (t r = 18) is similar to the low-amplitude case. But as the structure 
rotates its major axis toward the horizontal, a tail of weak vorticity is formed (see 
Figure 3 at t = 32). Successive snapshots show that the vortices continue to rotate, 
becoming vertical again when A\o reaches a second maximum. This is similar to 
the observations by MMB of active structures at higher Reynolds number. 

In Case 2D-3, we observe (see Figure 4) pairing of negative sign vorticity (in the 
lower portion of the layer) and shredding of positive sign vorticity (in the upper 
portion). The development time of the subharmonic is slow relative to the funda- 
mental, and pairing occurs around the time when the amplitude of the subharmonic 
reaches its maximum. 

8.2 Three-dimensional wake 

To study the sensitivity of the two-dimensional wake to initial (or inlet) three- 
dimensional disturbances, five cases were studied: two cases (Case 3D-1, and Case 
3D-1-1) where a two-dimensional mode plus background random noise of different 
amplitudes are used for initial conditions, and three cases with a two-dimensional 
mode and a pair of oblique modes (Case 3D-2-1, Case 3D-2-2 and Case 3D-2-3) of 
different amplitudes for initial disturbances. The relative phase between the two- 
dimensional mode and the oblique modes was kept the same at 7r/4 for all cases. 

The intent in Case 3D-1 is to mimic the experiment of Sato & Onda (1970) 
where a fundamental mode is forced, but, as in all experiments, background noise 
exists. Figure 5(a) shows the evolution of Aio for Case 2D-1 and 3D-1. Clearly 
the weak background noise in this case has little effect on the evolution of the two- 
dimensional layer. The structure that dominates in this case is the two-dimensional 
double roller. In this case, the fastest growing three-dimensional modes are those 
whose wave-numbers are closest to that of the two-dimensional fundamental. Large 
amplitude disturbances were used in Case 3D-1-1 (A\q — 0.1, average amplitude of 
three-dimensional modes of 0.025). In this case, a band of wave numbers around 
the fundamental are enhanced after the roll-up. But, the energy in the subharmonic 
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Figure 3. Spanwise vorticity in Case 2D-2 at t = 32; contour increments are 0.09 
and tic marks are at bij 2 intervals. 
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Figure 4. Spanwise vorticity in Case 2D-3 during pairing/ tearing. Contour incre- 
ments are 0.05 and tic marks are at 2 x b\t 2 intervals. 
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FIGURE 5. Evolution of the two- dimensioned fundamental mode amplitude in 
3D-1 and 2D-1 (plots on top of 3D-1 case). 


mode decreases after the rollup, suggesting that pairing in the wake is not as natural 
as in the mixing layer. But as seen in our two-dimensional simulations, it will occur 
when forced. The vortex structure in case 3D-1-1 is again dominated by the two- 
dimensional vortices. 

The results of the previous two cases suggest that low- wavenumber oblique modes 
are not amplified during the rollup of the fundamental structure in the wake. We 
thus use a relatively small spanwise wavelength ( X z = 0.6A X ) for the forced oblique 
modes as in Rogers &: Moser (1992). The three oblique cases; Case 3D-2-1, 3D 2-2 
and 3D- 2-3 have initial amplitudes such that A\ o > A\\, A\$ < A\\ , and Aio = A\\ 
respectively. 

In Case 3D-2-1, the evolutions of the mode amplitudes indicate that the two- 
dimensional mode grows, saturates, and decreases while the energy in the oblique 
modes reaches a first peak at the same time as the fundamental (tf = 14) but later 
continues to increase until a second maximum is reached at t = 30. At this time, the 
oblique mode amplitude is at almost the same level as the fundamental. Constant 
spanwise vorticity surfaces for this case at t = 14 and t = 30 are shown in Figure 6. 
In Figure 6(a), the core rolls are apparent but are surrounded by ‘wisps’ of vorticity 
that form hoop-like structures. Later (t = 30) after the oblique modes grow (at 
t = 30) we find that the strong spanwise vorticity is concentrated in compact regions 
(Figure 6(b)). At this time, the streamwise and the cross-stream components of the 
vorticity dominate the large scale structures. The counter-rotating vortices shown 
in Figure 6(c) become the dominant feature of the flow. Of the three vorticity 
components, u> y is largest in magnitude, while w z is smallest. The vortex structures 
are inclined at 40 to 50 degrees with respect to the streamwise direction. These 
streamwise vortex pairs are reminiscent of the rib vortices in the plane mixing 
layer. 
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FIGURE 6(a). Three-dimensional structure of spanwise vorticity in Case 3D-2-1 
at t — 14: Surfaces represent 50% of peak value (±2.6), and tic marks are b x / 2 
intervals. 



FIGURE 6(b). Three-dimensional structure of spanwise vorticity in Case 3D-2-1 
at t = 30: Surfaces represent 50% of peak value (±1.6), and tic marks are b X j 2 
intervals. 


When the oblique modes are initially dominant (Case 3D-2-2), it was found that 
they grow rapidly and suppress the growth of the two-dimensional fundamental 
mode. Since the oblique modes are larger than the two-dimensional mode, the 
spanwise vorticity structure is not observed in the wake. When all the modes have 
the same amplitude (Case 3D-2-3), the energy in the modes indicates that the 
oblique modes and the fundamental mode grow together and that the oblique mode 
saturates first, then the fundamental modes saturate. 


332 


H . Maekawa, R. D, Moser & N. N. Mansour 



FIGURE 6(c). Three-dimensional structure of streamwise vorticity in Case 3D-2-1 
at t = 30: Surfaces represent 50% of peak value (±1.9), and tic marks are b x / 2 
intervals. 



FIGURE 6(d). Three-dimensional structure of in Case 3D-2-1 at t = 30: Surfaces 
represent 50% of peak value (±2.6), and tic marks are 6 j /2 intervals. 


4* Conclusions 

* Direct numerical simulations have been used to study the three-dimensional evo- 
lution of a plane wake. The flow has been analyzed by visualizing the vortical 
structures using vorticity contours and by tracking the modal energy in time. The 
two-dimensional features are similar to the results obtained for a two-dimensional 
spatially developing wake (MMB). 
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The growth of fundamental mode generates characteristic alternating sign vor- 
tices in the wake. When a two-dimensional fundamental mode is forced with high 
amplitude, vortical structures with ‘wisps’ around the roller are generated. When 
the wake is forced with a fundamental and a subharmonic asymmetric pairing and 
shredding are observed. 

When weak three-dimensional random disturbances are imposed in addition to a 
two-dimensional fundamental mode, we find that the two-dimensional rollers appear 
to be the dominant features of the wake. Analysis of the energy contained in each 
mode reveals that three-dimensional modes with wavenumber close to that of the 
two-dimensional fundamental mode grow after saturation of the fundamental mode. 
In contrast, the energy in the low wavenumbers modes decreases. 

When the wake is forced with oblique modes of various amplitudes, we find that 
at large amplitudes the growth of the oblique mode leads to the breakdown of 
the alternate-sign vortex structures in the wake, and counter-rotating streamwise 
vortical structures are generated after the breakdown. These streamwise structures 
are inclined at 40 to 50 degrees with respect to the streamwise direction. When the 
amplitude of the initial oblique modes is higher than that of the fundamental, the 
growth of the fundamental mode is suppressed, and no spanwise coherent structures 
were found in this case. 
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